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Abstract. We present an alternative method for the kinematic analysis of high proper motion surveys and discuss 
its application to the survey of Oppenheimer et al. (2001) for the selection of reliable halo white dwarfs (WDs). 
The local WD space density we estimate is pwD — 1 -j- 2 • 10 -5 MqPC -3 , which is about an order of magnitude 
smaller than the value derived in Oppenheimer et al. (2001), and is consistent with the values obtained from 
recent reanalyses of the same data (e.g. Reid et al. 2001, Reyle et al. 2001, Torres et al. 2002, Salim et al. 2004). 
Our result, which corresponds to a fraction of 0.1% -5- 0.2% of the local dark matter, does not support the scenario 
suggested by the microlensing experiments that ancient cool WDs could contribute significantly to the dark halo 
of the Milky Way. 
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1. Introduction 

One of the most recent challenges in observational astron- 
omy is to explain the nature of the objects that produced 
the microlensing events towards the Magellanic Clouds 
(Alcock et al. 2000). The most obvious candidates for 
these events are ancient white dwarfs, so that several 
projects have been carried out in recent years to reveal 
the existence of such hidden population of dim sources 
(see Hansen & Liebert 2003 for a review). The most ex- 
tensive survey to date is that of Oppenheimer et al. (2001, 
OHDHS). They discovered 38 suspected halo white dwarfs 
and derived a local density of ~ 1.1 ■ 10~ 4 M Q pc -3 , which 
corresponds to a fraction of 1-2% of the halo dark mat- 
ter in the vicinity of the Sun. Different authors challenged 
these results on the basis of the age estimates of the candi- 
dates (Hansen 2001, Bergeron 2003), or after a reanalysis 
of the kinematic data (e.g. Reid, Sahu & Hawley 2001; 
Reyle et al. 2001, Flynn et al. 2003, Torres et al. 2002). 

In any event, all those studies evidence a significant 
contamination of thick disk objects affecting the halo WD 
sample, and point out the basic problem of defining an 
accurate procedure to deconvolve the halo and thick disk 
populations on the basis of their kinematic and photomet- 
ric properties. 

In this paper we describe a general statistical method 
designed to reject objects with disk kinematics and iso- 
late probable halo members from the screening of kine- 
matically selected samples. Finally, we discuss the results 



obtained with this method when applied to the OHDHS 
survey, and compare them to the preliminary results de- 
rived from the GSC Il-based new high proper motion sur- 
vey in the Northern hemisphere by Carollo et al. (2004). 



2. SSS Halo WD survey 

The OHDHS survey was based on digitized, photo- 
graphic Schmidt plates (R59F and Bj passbands) from 
the SuperCOSMOS Sky Survey (SSS, Hambly et al. 2001). 
They analyzed 196 three epoch plates (IllaJ, IIIaF and 
IV-N) covering an area of 4165 square degrees near the 
South Galactic Pole (SGP). The magnitude limit of the 
survey is of R59F = 19.8, while the proper motion lim- 
its are 0.33"yr -1 < H < 10"yr -1 - They found 98 WDs, 
whose tangential velocities were derived from the mea- 
sured proper motions and photometric distances esti- 
mated via a linear color magnitude (CM) relation, Mb } 
vs. Bj — R59F, calibrated by means of the WD sam- 
ple with available trigonometric parallaxes published by 
Bergeron, Ruiz & Legget (1997). The kinematic analy- 
sis of this sample was made in the two dimensional (U,V) 
plane, after assuming that the third galactic velocity com- 
ponent was zero (W — 0). Thick disk contaminants were 
rejected with a 2cr threshold, ^U 2 + (V + 35) 2 > 95 km 
s _1 , which would correspond to a 86% confidence level 
in the case of a non-kinematically selected sample. In 
this way, 38 WDs were considered as halo members, from 
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which a space density of pwd — 1.1 1CP 4 M Q pc~ 3 was 
computed, assuming 0.6 M Q for the average WD mass. 

As mentioned in the previous section, these results 
were critically revised by several authors. In particular, 
an independent kinematic analysis of the OHDHS sam- 
ple was performed by Reid et al. (2001), who noted that 
the resulting distribution of the WDs in the (U,V) dia- 
gram seems more compatible with the high velocity tails 
of the thick disk. They computed (U, V) components as- 
suming that the unknown radial velocity is null (V r — 0) 
and selected halo WDs with the crude but robust crite- 
rion of accepting objects with retrograde motion only (4 
objects). This leads to a more conservative value of the 
density, p WD ~ 1.8 10~ 5 M pc~ 3 . 

Recently, Salim ct al. (2004) reanalyzed the WD sam- 
ple of OHDHS on the basis of new spectroscopic and pho- 
tometric measurements. Radial velocities of 13 WDs with 
H a line, and standard Johnson-Cousins photometry for 
half of the sample were obtained. In addition, distances 
were redetermined with the CCD photometry by means 
of the theoretical color magnitude relation for hydrogen 
and helium atmospheres published by Bergeron, Leggett 
& Ruiz (2001). Salim et al. (2003) confirmed the results 
of OHDHS with the same 95 km s" 1 (2a) threshold, but 
showed that a minimum density, nwD — 3.1 10 -5 pc~ 3 
is attained with a higher, more conservative, threshold of 
190 km s _x . 



3. Kinematic analysis 

The kinematic analysis of the WD sample drawn from 
a proper motion limited survey, including the choice of 
a criterium for rejecting the contaminant disk WDs and 
select the true halo WDs, is one of the critical steps of this 
kind of studies. 

As the velocity distribution of the disk(s) and halo 
population do partially overlap (Fig. 1), it is not possible 
to infer univocally, on the basis of kinematic data alone, 
the parent population of every object. Nevertheless, it is 
always possible to test if an object is, or is not, consistent 
with the velocity distribution of a certain population once 
a value for the confidence level is chosen. 

Here, we retain as halo WDs those objects whose kine- 
matics is not consistent with the velocity distribution of 
the thick disk population 1 given a certain confidence level; 
this allows the identification of halo WDs while limit- 
ing the contamination of high velocity thick disk objects. 
Unless corrected for the incompleteness due to the frac- 
tion of rejected halo WDs whose kinematics is compatible 
with that of the thick disk population, it is clear that this 
procedure can only provide a lower limit to the actual 
density. 

An alternative, and potentially more rigorous proce- 
dure, is a Maximum-likelihood analysis that fits simulta- 
neously the superposition of two or more populations (see 



1 Implicitly, we assume that besides the thick disk WDs, this 
criterion rejects the "slowest" thin disk objects as well. 



e.g. Nelson et al. 2002, Koopmans & Blandford 2002). In 
this case however, because of the small size of the samples, 
further assumptions on the kinematics and the formation 
process (IMF, age, etc.) of all the populations involved are 
usually necessary. 

3.1. Schwarzschild distribution 

We assume that the probability that the galactic velocity 
components (U,V,W) of an object in the solar neighbor- 
hood belonging to a certain stellar population lies in the 
element of velocity space d 3 v — dU dVdW is well described 
by a Schwarzschild distribution: 
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which represents a trivariate gaussian ellipsoid, where 
Vo indicates the rotation lag with respect to the LSR and 
cry, ay, and aw the velocity dispersions. 

In practice, the galactic components need to be derived 
from the observed tangential and radial velocity compo- 
nents (V Q , Vg, V r ): 



(2) 



where G2000 = G(a,S) is the transformation matrix 
from the equatorial coordinates system (J2000) to the 
galactic system, which depends explicitly on the stellar 
position (a, S). Here, (Uq,Vq,Wq) is the Sun veloc- 
ity with respect to the Local Standard of Rest (LSR), 
for which Dehnen & Binney (1998) estimated (+10.00 ± 
0.36, +5.25 ± 0.62, +7.17 ± 0.38) km s" 1 from the analy- 
sis of the Hipparcos catalogue. The tangential velocities 
V a and Vg (km s _1 ), are computed from the observed 
proper motions (arcsec yr _1 ) and distances (pc) derived 
from trigonometric or photometric parallaxes, ir = l/d, as 
usual: 

V a = 4.74047 dp a cos 5 
V 5 = 4.74047 d p 5 

3.2. Tangential velocity distribution 

If the full 3D space velocity cannot be recovered, as in 
the case of proper motion surveys, we can adopt a similar 
procedure in the 2D tangential velocity plane, (V a , Vs). 
The bivariate marginal distribution, ip{V a , Vs), can be ob- 
tained by properly integrating the distribution in Eq. ^ 
along the V r component: 
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Fig. 1. Tangential velocity distributions, ip(V a ,Vs), to- 
ward the direction (/ = 143°, b = -57°) of WD0135- 
039 (solid circle). The ellipses show the iso-probability 
contours (2a, 3a) of the thin disk, thick disk, and lo- 
cal halo populations based on the kinematic parameters 
from Binney & Merrifield (1998), Soubiran et al. (2003) 
and Casertano, Ratnatunga & Bahcall (1990), respec- 
tively. The three concentric circles indicate the velocity 
thresholds, V m i n = 4.74^iii m r, for the OHDHS survey 
(Mlim = 0.33" yr _1 ) at the distances r = 40 pc, 80 pc 
and 120 pc. 



This is a general bivariate gaussian distribution which is 
defined by five parameters: V a o, Vao, ay a, &vs and p. 
These parameters are linear functions of the first and sec- 
ond order moments of Eq. ^ as described for instance in 
Trumpler & Weaver (1953). 

Our analysis will be based on Eq. 3 that represents the 
appropriate density distribution when radial velocities are 
missing. 

Notice that this approach, even in the case of surveys in- 
volving widely different line-of-sights, allows the deriva- 
tion of the exact tangential velocity distribution for every 
star, without any assumption on the unknown third ve- 
locity component V r . 

3.3. Thick disk model 

The following properties for the population of thick disk 
WDs in the solar neighborhood were assumed: 

— a uniform local space density; for, the typical distance 
reachable by ground based surveys (~ 100 pc) is much 
smaller than the exponential vertical scale-height of 
the thick disk (h z ~ 1000 pc); 



— a velocity distribution (Eq.[TJ with (ajj, ay, aw, Vq) ~ 
(63,39,39,-45) km s _1 , as derived in Soubiran, 
Bienayme & Siebert (2003). 

We notice that the velocity ellipsoid of the thick disk 
population is not currently well established so that this 
choice will somehow affect the final result. For instance, 
the presence of a non-gaussian high velocity tail (cfr. 
Gilmore et al. 2002) would increase the contamination af- 
fecting the halo WD sample. 

3.4. Kinematically selected samples 

In the case of a magnitude- and /x-limited survey with a 
total extension of £1 steradians, the following observational 
constraints need to be taken into account: 

1. an apparent magnitude limit m < mn m which implies a 
distance limit as a function of the absolute magnitude, 
M, of the target: 

r < W(M) = l [0-2( m n m -M) + l] 

which, in turn, defines the maximum volume 2 , 
V Max (M) = |r2rf nax covered by the survey; 

2. a proper motion limit /i > /xu m which translates into a 
tangential velocity threshold varying with stellar dis- 
tance: 

V tan = Jv[+vj > V mia (r) = 4.74 Wlm r. 

Note that, although the distance distribution (oc r 2 ) 
and the kinematic distribution (Eq. 3 ) of the complete 
population are independent, now they result correlated 
for the observed sample because of the existence of the 
velocity threshold, V m i n (r). 

The probability to select a star with absolute magni- 
tude M in the range (r, r + dr), (V a , V a + dV a ), (Vs,Vs + 
dV$) is then dP = f(r,V a ,Vs)drdV a dVs, where the joint 
probability density is: 

f Kr 2 ^(V a ,V 5 ) if V tan > V min (r) 

n t/ t/ \ J and r<r max (M) 

/(r '^ )= ifl4 an <Wr) W 

[ or r > r max (M) 

Here, K is a normalization constant such that 
fj J fdrdV a dV 5 - 1. 

If we integrate over r the joint probability density func- 
tion given in Eq. we obtain the marginal density distri- 
bution 

h(V a ,V s ) = I f(r,V a ,V s )dr (5) 
Jo 

2 Note that this is a purely photometric definition which 
does not correspond exactly to the analogue quantity adopted 
for the evaluation of the WD density via the 1/V ma x method 
(Schmidt 1975). 
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which quantifies the probability that an object with tan- 
gential velocities (V a , Vg) could be randomly found some- 
where within the whole volume if2rj^ ax , where an ob- 
ject with absolute magnitude M could in principle be ob- 
served. 

At the same time, we can introduce the (condi- 
tional) probability that an object with tangential veloc- 
ities (V a , Vg) is found at the measured distance, r: 

,/ T/ t 7 1 \ f(r,V a ,V s ) 

t(V a ,Vs\r) = — (6) 

9(r) 

where 

9{r)= / / f(r,V a ,V s )dV a dV 5 (7) 

J J — oo 

is the marginal density distribution which defines the 
probability that an object with whatever velocity can be 
observed at a distance r. Because the velocity threshold 
increases linearly with distance, Kni n oc r, the space dis- 
tribution of the proper-motion selected sample (Eq. 7) is 
also biased towards smaller distances. 

Both the marginal distribution h(V a , Vg) and the con- 
ditional probability t(V a , Vg \r) can be used to test the con- 
sistency of each object with a parent population. In prin- 
ciple, the conditional probability t(V a ,Vi|r) seems more 
appropriate than h(V a ,V$) since it fully utilizes the indi- 
vidual stellar distances. However, the differences become 
insignificant when the confidence level is set to sufficiently 
high values (see next section). 

Note that, formally, Eq. El is equivalent to the original 
distribution, ip(Va, Vg), except that the probability is null 
for y/V£ + Vg < V ni i n (r), and it has been re-normalized. 

3.5. Confidence intervals and contamination 

Basically, because a proper motion limited survey under- 
samples the low velocity objects, the main difference be- 
tween the kinematically selected distributions (Eqs. 15161) 
and the complete one (Eq. 3) is that the probability den- 
sity is redistributed from the low velocity regions towards 
the high velocity tails. This means that the observed sam- 
ple is biased towards high velocity objects, as shown for in- 
stance by the simulations of Reyle et al. (2001) and Torres 
et al. (2002). 

This effect needs to be taken into account when we 
define a confidence interval over the (V a ,Vg) plane in or- 
der to test the consistency with the parent population and 
to estimate the contamination due to objects in the tails 
beyond the critical limit. In fact, the adoption of the orig- 
inal ip(V a ,Vg) to reject the disk stars with respect to a 
certain confidence level, e.g. 1 — a = 99%, would exclude 
99% of all the existing thick disk stars which, however, 
corresponds to a smaller fraction of the thick disk objects 
that are really present in the kinematically selected sub- 
sample. In this case, only the confidence interval defined 
for t(V a ,Vg\r), or h(V a ,Vg), assures that the fraction of 
false negatives contaminating the sample of bona fide halo 



stars does not exceed - on average - 1% of the observed 
thick disk objects. 

In the left panels of Figures 12131 the concentric el- 
lipses show the iso-probability contours (ler, 2a, 3a) of 
the velocity distribution expected for thick disk stars, 
ip(V a , Vg), evaluated in the direction of one of the stars in 
the Oppenheimer's sample (LHS 1447), whose tangential 
velocity is marked with a filled circle. The points represent 
a Montecarlo realization of 2000 simulated WD's drawn 
from the kinematically selected distributions h(V a ,Vg) 
and t(V a , Vg\r). The excess of "simulated" thick disk stars 
with high velocity is evidenced by the fact that there are 
many more than ~20 objects (1% of the simulated sample) 
outside the 3a confidence interval. 

LHS 1447 is also located outside the 3a contour so 
that, according to the complete distribution, it should be 
rejected as a thick disk star with a confidence level higher 
than 1— a =99%. Actually, that conclusion would be incor- 
rect if we tested the hypothesis that LHS 1447 is a member 
of the kinematically selected sample as shown in the right 
panels of Figures 12131 where the marginal and conditional 
distributions, h(V a , Vg) and t(V a , Vg\r), are drawn. In fact, 
in these cases the star is located within the iso-probability 
contour delimiting the 99% confidence level so that it must 
be accepted as a thick disk star. 

4. Results and discussion 

Both distributions, h(V a , Vg) and t(V a ,Vg\r), were used to 
analyze the WD sample in the OHDHS survey. 

The kinematic tests were carried out in the tangential 
plane of each individual star so that no assumption on ra- 
dial velocity is necessary. The values of 95% and 99% for 
the confidence level (1 — a) were chosen in order to mini- 
mize the presence of false negatives. With a total sample 
of 98 WDs, presumably a mixture of (thin and thick) disk 
and halo WDs, we expect that < 1 (99%) and < 5 (95%) 
of the high velocity thick disk stars would contaminate the 
selected Pop. II WD sample. 

4.1. Halo WD density 

In Table Q we report the results based on this proce- 
dure for the WD sample published by OHDHS. We only 
found 10 objects which do not appear consistent with the 
kinematically selected density distributions, h(V a , Vg) and 
t(V a , Vg\r), at the 99% confidence level, while 12-13 proba- 
ble halo WDs are selected when l — a = 95%. As expected, 
the number of candidates increases up to 14 (99%) or 20 
(95%) in the case of a test based on the complete distribu- 
tion, tp(V a , Vg), mainly because of a higher contamination. 

Finally, the halo WD density was estimated by means 
of the classical l/Vjyiax method (Schmidt 1975), and as- 
suming a value of 0.6 M Q for the typical WD mass. The 
results, with their (poissonian only) errors, are reported 
in Tab. ^ where the different values refer to the two con- 
fidence levels and the three probability distributions used 
for the calculations. 
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Fig. 2. Left panel: iso-probability contours (la, 2a, 3a) oiip(V a , Vs) compared against a Montecarlo simulation (dots) 
of the thick disk stars in the direction (I = 226?34,6 = — 64?27) of LHS 1447 (solid circle) drawn from the kinematically 
selected distribution, h(V a ,Vs). Right panel: iso-probability contours (confidence levels of 40%, 70%, 95% and 99%) 

of h(y a , v s ). 
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Fig. 3. Left panel: iso-probability contours (la, 2a, 3a) oiip(V a , Vs) compared against a Montecarlo simulation (dots) 
of the thick disk stars in the direction (I = 226?34,6 = — 64?27) of LHS 1447 (solid circle) drawn from the kinematically 
selected distribution, t(V a , Vs\r). Right panel: iso-probability contours (confidence levels of 40%, 70%, 95% and 99%) 
oft(V a ,V*|r). 



Although affected by large uncertainties, the values in 
Tab. 2] suggest a density of p W D ~ 10~ 5 M pc -3 , i.e. 
0.1-0.2% of the local dark matter, which is an order of 
magnitude smaller than what reported in OHDHS. 

Our results are consistent with the local mass den- 
sity of halo WDs estimated by Gould et al. (1998), and 



with various reanalyses of the OHDHS sample (e.g. Rcid 
et al. 2001, Reyle et al. 2001, Torres et al. 2002, Salim 
et al. 2004). Furthermore, Carollo et al. (2004), applying 
the statistical methodology described in this paper on a 
new high proper motion survey based on GSC-II material, 
derived a similar value of ~ 10~ 5 M Q pc~ 3 . 
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Table 1. Estimation of the halo WD density based on the objects selected from the OHDHS sample. 



Confid. 




i>{v a ,v s ) 




h(V a ,V s ) 




t(V a ,Vs\r) 


level 


WDs 


Pwd (M pc- 3 ) 


WDs 


Pwd (M pc- 3 ) 


WDs 


Pwd (M pc" 3 ) 


99% 


14 


(2.0 ± 0.9) • 10 _b 

/n i I i o\ i n-5 


10 


(1.6 ±0.8) • 10 _b 


10 


(1.5 ±0.8) • 10 _b 



95% 20 (3.1 ± 1.0) ■ 10~ b 12 (1.9 ±0.9) ■ 10~ b 13 (2.0 ± 0.9) ■ 10~ b 



Table 2. Estimation of the halo WD density based on the sample revised by Salim et al. (2004). 



Confid. 








h(V a ,Vs) 




t(V a ,Vs\r) 


level 


WDs 


Pwd (M pc 3 ) 


WDs 


Pwd (M pc- 3 ) 


WDs 


Pwd (M Q pc -3 ) 


99% 


19 


(1.8 ±0.6) • 10" b 


17 


(1.6 ±0.6) • 10"" 


16 


(1.5 ±0.6) • 10"" 


95% 


28 


(3.3 ± 0.8) • 10" b 


18 


(1.8 ±0.6) • 10~ b 


18 


(1.9 ±0.6) • 10~ b 



Lacking individual trigonometric parallaxes, a criti- 
cal point of this (and any) analysis is the choice of the 
method for the estimation of the distances, which directly 
affects the evaluation of the WD tangential velocities and, 
of course, of their stellar density. As remarked by several 
authors (see e.g. Torres et al. 2002, Bergeron 2003), em- 
pirical and theoretical CM relations can both give rise to 
systematic errors. 

To this regard, if for the distances of the OHDHS sam- 
ple we adopt the values recently redetermined 3 by Salim 
et al. (2004), the number of selected halo WDs increases 
but the resulting densities, shown in Table [21 are not sig- 
nificantly different from those reported in Table 

4.2. Distance and velocity errors 

The large error, ~ 20-30%, affecting WD photometric par- 
allaxes, cannot be neglected in a rigorous statistical anal- 
ysis. Basically, besides the contribution of the photomet- 
ric errors, the large uncertainty in the distance modu- 
lus, m — M, derives from the large intrinsic dispersion 
(<7Mv — 0.4 - 0.5 mag) of the CM relation, a consequence 
of the superposition of cooling sequences of WDs of dif- 
ferent masses and atmospheres. 

In practice, the main effect of the tangential velocity 
errors, ey/V = \J (cr^/fi) 2 + (<Td/d) 2 , is to increase the 
dispersion and the overlap of the "observed" kinematic 
distributions belonging to the various stellar populations. 
Clearly this also increases the contamination of the disk 
WDs and makes the identification of the halo WDs more 
difficult. 

Although a more rigorous statistical analysis should 
be necessary to consider properly the presence of these 
errors, a conservative estimation can be given by select- 
ing only those objects which are not consistent with the 
"observed" kinematic distribution that results from con- 
volving the projected kinematic distribution of the thick 
disk (Eq. 3) with a bivariate gaussian error distribution 
with null mean and dispersions, (ey Q , evs)(i), correspond- 

3 They adopted CM relations based on theoretical cooling 
tracks of 0.6 M WDs with H or He atmospheres. This resulted 
in distances 16% systematically larger (on average) than those 
in OHDHS. 



ing to the velocity errors of the i-th object. The velocity 
errors have been derived by assuming the proper motion 
errors, er^, listed in Tab. 1 of OHDHS, and a more real- 
istic photometric parallax error, ad/d, of 25% (instead of 
20%). 

The different halo WD densities estimated from 
the objects which are not consistent (at the 95% and 
99% confidence level) with the new distributions are 
reported in Tables 3 and 4. Because of the larger velocity 
thresholds, the number of selected halo WDs is smaller 
than those reported in Tables 1 and 2. The estimated 
WD densities, uncorrected for the loss of halo WDs with 
disk kinematics, decrease proportionally, but are still 
consistent with pwd ~ 10 -5 M Q pc -3 . Note that the 
minimum values, which are reported in Tab. 4, have been 
derived from the data of Salim et al. (2004) who provided 
distances (and thus volumes) systematically larger than 
OHDHS. 



4.3. On the thick disk model 

As mentioned in Sect. 13.31 our selection criterion depends 
implicitly also on the choice of the kinematic parameters 
adopted for the thick disk, whose spatial and kinemat- 
ical properties are still matter of debate and investiga- 
tion. Here, we have used the velocity ellipsoid recently 
derived by Soubiran et al. (2003) from a sample of ~ 400 
giants with 3D kinematics at a distance of 200-800 pc to- 
wards the North Galactic Cap. Their results are very close 
to the kinematic parameters estimated 4 by Casertano, 
Ratnatunga & Bahcall (1990) and are consistent with var- 
ious other determinations of the thick disk kinematics, 
which support velocity dispersions of 40-60 km s^ 1 and 
an asymmetric drift in the range 30-50 km s _1 . 

Although controversial, some authors claim the pres- 
ence of a vertical velocity gradient, that supports a thick 
disk which rotates faster close to the galactic plane (i.e. 
where the WD sample is localized), than at higher Z's, 

4 Casertano, Ratnatunga & Bahcall (1990) derived 
(<tu,<tv,<tw,Vo) ^ (66,37,38,-40) ± 10 km s" 1 from a 
maximum likelihood analysis of high proper motion stars 
within 500 pc of the Sun. 
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Confid. ip(V a ,V s ) h{V a ,V s ) t(V a ,V s \r) 

level WDs p W p (M Q pc~ 3 ) WDs p W p (M pc~ 3 ) WDs p W p (M pc~ 3 ) 

99% 6 (1.3 ±0.8) ■ l(T b 3 (1.2 ±0.8) ■ l(T b 3 (1.2 ± 0.8) ■ 10~ b 

95% 14 (2.0 ±0.9) ■ 10~ b 5 (1.3 ± 0.8) ■ 10~ b 6 (1.3 ± 0.8) ■ 10~ b 



Table 4. Same as Tab|5] after adopting a thick disk velocity distribution convolved with the observation errors. 



Confid. i>{V a ,Vs) h(V a ,Vs) t{V a ,Vs\r) 

level WDs p W p (M pc~ 3 ) WDs p W p (M pc~ 3 ) WDs p W p (M pc~ 3 ) 

99% 14 (1.2 ±0.5) ■ 10~ b 8 (0.7 ± 0.4) ■ 10~ b 5 (0.6 ± 0.4) ■ 10~ b 

95% 18 (1.7 ±0.6) ■ 10~ b 14 (1.2 ± 0.5) ■ 10~ b 9 (1.0 ± 0.5) ■ 10~ b 



where the studies of the thick disk kinematics have been 
usually carried out. In particular, Chiba & Beers (2000), 
who analyzed 1203 metal poor stars non-kinematically se- 
lected, found a rapidly rotating thick disk close to the 
galactic plane with a small asymmetric drift Vq ~ —20 
km s _1 and with velocity dispersions {ou, °~v , °~w) — 
(46 ± 4, 50 ± 4,35 ± 3) km s -1 . Moreover, they deter- 
mined a velocity gradient dVo/d\Z\ ~ —30 ± 3 km s _1 
kpc~\ that, however, other studies (e.g. Soubiran et al. 
2003) do not detect. Nevertheless, a fast rotating thick 
disk at Z was determined 5 also by Upgren et al. 
(1997) from a sample of K-M dwarfs in the solar neigh- 
borhood (d <> 50 pc) with trigonometric parallaxes and 
proper motions from the Hipparcos catalogue and radial 
velocity measurements. 

Thus, in order to test the sensitivity of our method 
with respect to the adopted thick disk model, we repeated 
the WD selection of the Salim et al. (2003) sample through 
the distributions h(V a ,Vs) and t(V a ,Vs\r) derived using 
the velocity ellipsoid from Chiba & Beers (2000). The new 
results are consistent (within la) with the values obtained 
with the kinematics from Soubiran et al. (2003), although 
the resulting densities appear typically larger than the pre- 
vious ones. 

For instance, with a 99% confidence level we find 
p WD ^ (1.7 ± 0.6)10" 5 M pc~ 3 for both h(V a ,Vg) and 
t{V a ,Vs\r) when the velocity errors are not taken into 
account (cfr. Tab. |2J, while the distributions convolved 
with the velocity errors provide pwd — (0.9 ± 0.5)10 -5 
M Q pc~ 3 (cfr. Tab.E|). The 95% confidence level also pro- 
vides similar but systematically higher new densities up 
to (1.8 ±0.7)10" 5 M pc- 3 and (3.0 ± 0.9)10~ 5 M Q p C - 3 
respectively when the velocity errors are, or are not, con- 
volved with the tangential velocity distributions. 

Anyhow, it appears that, with the adopted confidence 
levels, significantly higher density (e.g. close to ~ 10~ 4 
Mqpc -3 may be attained only with disk ellipsoids kine- 
matically much "cooler" than those expected for a typ- 
ical thick disk population. For instance, a total density 
(8.8 ± 0.2) 10~ 5 M pc~ 3 is only obtained counting all the 

5 They estimated a rotation lag of Vb — —28.3 ± 3.8 km s _1 
for the "old" disk component with dispersions (cry, ov, <?w) — 
(56.1 ± 3.9, 34.2 ± 2.5, 31.2 ± 2.5) km s" 1 . 



41 WDs which are not consistent with the thin disk 6 kine- 
matics (using t(V a , Vs\r) with a 95% confidence level), i.e. 
summing both halo and thick disk WDs. 

4.4. UVW distribution 

Salim et al. (2003) provide radial velocities for 15 DA 
WDs, 13 of which derived from new measurements of the 
OHDHS sample and two from Pauli et al. (2003), so that, 
in principle, a more accurate kinematic membership for 
these objects may be inferred using the information from 
the full 3D velocities. This requires 3D velocity distribu- 
tions for kinematically selected samples which are beyond 
the scope of the current study. However, the availability 
of both tangential and radial velocities for this subsample 
offers the possibility to check a posteriori the efficiency 
of the 2D kinematic analysis adopted in this work and 
described in Sect. [31 

To this regard, Figure [5] shows the (U,V,W) velocities 
derived from Eq. |3 for the 15 stars with available radial 
velocity. Those which have been selected with a 95% con- 
fidence level by means of the distributions h(V a ,Vs) and 
t(V a ,Vs\r) convolved with the velocity errors (Tab.0J are 
marked with square and diamond symbols. In addition, 
the 3a iso-probability ellipses of the thick disk and halo 
velocity distributions, based on the kinematic parameters 
respectively from Soubiran et al. (2003) and Casertano, 
Ratnatunga & Bahcall (1990), are also plotted. The three 
panels of Fig. [3 indicate that, basically, all the likely halo 
WDs have been properly identified by our kinematic anal- 
ysis based on the 2D (V a , Vg) distributions, thus support- 
ing the reliability of our selection procedure. 

5. Conclusions 

A kinematically selected sample made of 98 WDs with 
// > 0.33"yr~ 1 was published by OHDHS who performed a 
high proper motion survey over 4165 deg 2 toward the SGP 
down to R59F ~ 19.8. These data stimulated a number 
of studies addressing the issue that a significant part of 
the dark halo of the Milky Way could be composed of 

6 We adopted (a n , ov, a w , V ) ^ (34, 21, 18; -6) km s" 1 
from Tab. 10.4 of Binney & Merrifield (1998). 
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Fig. 4. Velocity distribution (U,V,W) of the subsample of 15 stars with available radial velocity (dots with la error 
bars) from Salim et al. (2004). The objects selected by means of the distributions h(V a , Vg) and t(V a , Vs\r) with a 
95% confidence level (Tab. 0J) are marked with circle and diamond symbols, respectively. The ellipses show the 3cr 
iso-probability contours of the thick disk and halo velocity distribution. 



matter in the form of ancient cool WDs. The basic problem 
- as addressed by several authors - is the criterion to 
disentangle the mixture of (thick) disk and halo objects 
on the basis of their kinematic properties and ages. 

To this regard, we have implemented a general method 
for the kinematic analysis of high proper motion sur- 
veys and applied it to the identification of reliable halo 
stars. The kinematically-selected tangential velocity dis- 
tributions are derived for every star, so that no assump- 
tion on the unknown third velocity component, V r , nor 
any approximation on the galactic components (U,V,W), 
is necessary. 

We selected as bona fide halo WDs only those stars 
whose tangential velocity is inconsistent, at the 95% and 
99% confidence levels, with the appropriate projected dis- 
tribution, h(V a ,V$) or t(V a ,Vs\r), of the observed thick 
disk population, thus assuring limited contamination of 
thick disk objects. Finally, the effect of large velocity er- 
rors, which derive from the intrinsic uncertainty of the 
WD photometric parallaxes, was also discussed and taken 
into account. 

We applied this methodology to the OHDHS sample 
and selected 10 probable halo WDs (that became 3 after 
the inclusion of the velocity errors) at the a 99% confidence 
level. Through the 1/VMax method, we estimated a local 
WD density of p W D ^ 1-^2- 10~ 5 M pc- 3 (i.e. 0.1-0.2% of 
the local dark matter) which is consistent with the values 
found by Gould et al. (1998), as well as by other authors 
who reanalyzed the OHDHS sample (e.g. Reid et al. 2001, 
Reyle et al. 2001, Torres et al. 2002, Flynn et al. 2003). 
The same methodology applied to the OHDHS sample 
revised by Salim et al. (2004) yields a similar value. These 
results agree with those found by Carollo et al. (2004) 
from a first analysis of new data of an independent high 
proper motion survey in the Northern hemisphere based 



on material and procedures used for the construction of 
the GSC-II. 

Although affected by a large uncertainty due to the 
small statistics and low accuracy of the photometric par- 
allaxes, our results clearly indicate that ancient cool WDs 
do not contribute significantly to the baryonic fraction of 
the galactic dark halo, as possibly suggested by the mi- 
crolensing experiments which claimed that ~ 20% of the 
dark matter is formed by compact objects of ~ 0.5 M Q 
(Alcock et al. 2000). 
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